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, We study the phase diagram at finite temperature of Josephson junction 

arrays with capacitive disorder (i.e., random offset charges and/or random 
c/3 ' charging energies): in the limit of large particle numbers per junction, this 

is a remarkable physical realization of the disordered boson Hubbard model. 
By using a mean-field approximation, we compute the average free energy 
and the equation for the phase boundary line between the insulating and 
^ , the superconducting phase. We find that the Mott-insulating lobe structure 

O \ disappears for large variance {a ^ e) of the offset charges probability distribu- 

tion. Further, with nearest-neighbor interactions, the insulating lobe around 

■ q = eis destroyed even for small values of o". In the case of random charging 
^ , energies, until the variance of the distribution reaches some critical value the 

^ \ superconducting phase increases in comparison to the situation in which all 

\ self-capacitances are equal. 
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I. INTRODUCTION 



In practical realizations of Josephson devices [1], one has to deal with capacitance disorder 
caused either by offset charge defects in the junctions or in the substrate (random offset 
charges) [2] or by imperfections in the construction of the devices, which may lead to random 
capacitances of the Josephson junctions. Random offset charges cannot be made to vanish by 
using a gate for each superconducting island since in large arrays too many electrodes would 
be necessary (i.e., too complicated fabrication procedures). In principle, an external uniform 
charge can be introduced and tuned in Josephson junctions arrays (JJA's) by applying 
a gate voltage with respect to the ground plane: in Ref. [3] this situation was analyzed 
experimentally by placing a gate underneath a Josephson array and it was observed a sensible 
variation (~ 40%) of the resistance between the unfrustratcd and the fully frustrated array. 
Although from a theoretical point of view charge and magnetic frustrations are dual to each 
other, experimentally it is possible to tune only the magnetic frustration in a controlled 
way. For this reason, it is widely believed that a challenging task for the theory is to 
develop reliable techniques to investigate the effects of random charge frustration on the 
phase structure of the arrays. 

In this paper, we shall address the problem of determining the finite temperature phase 
diagram of JJA's with capacitive disorder (i.e., with random offset charges and/or random 
self-capacitances). To derive the phase boundary between the insulating and the supercon- 
ducting phase, we shall use a mean-field (MF) theory approach in the path-integral approach 
for quantum JJA's with offset charges and general capacitance matrices [4,5]. We find that 
the charge disorder supports superconductivity and that the relative variations of the insu- 
lating and superconducting regions depend on the mean value q of the charge probability 
distribution: when g = 0, increasing the disorder leads to an enlargement of the supercon- 
ducting phase. If the charge disorder is sufficiently strong (a > e), the lobe structure [1] 
disappears: in other words, the phase boundary line (and the correlation functions) does 
not depend any longer on q. In the following, we shall provide a quantitative analysis of 
this phenomenon. We shall consider Gaussian, uniform (as in Ref. [6]) and 5- like distribu- 
tions. For T = our results agree with the results obtained in Refs. [6] and [7]. Also the 
randomness of the self-capacitances leads to remarkable effects: the superconducting phase 
increases with respect to the case where disorder is not present. We shall also compare the 
results of the functional MF approach with the molecular MF discussed in Appendix B. 

The Hamiltonian commonly used to describe the Cooper pair tunneling in supercon- 
ducting quantum networks defines the so called quantum phase model (QPM). In its most 
general form it is given by 



where Lpi is the phase of the superconducting order parameter at grain i. Its conjugate 
variable n\ ([(/9i,ni] = describes the number of Cooper pairs in the ith superconducting 
grain. The symbol {ij) indicates a sum over nearest-neighbor grains only. The second term 
in the Hamiltonian (1) describes the hopping of Cooper pairs between neighboring sites {Ej 
is the Josephson energy). The first term determines the electrostatic coupling between the 
Cooper pairs: Qi is the excess of charge due to Cooper pairs (Qi = 2en\) on site i and is 
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the capacitance matrix. An external gate voltage Vi gives the contribution to the energy via 
the offset charge qi = X^jCyVj. This external voltage can be cither applied to the ground 
plane or, more interestingly, it may be induced by charges trapped in the substrate. In the 
latter situation qi is a random variable: the effects of this randomness on the phase diagram 
are the main object of our investigation. We shall also treat explicitly the case of random 
self-capacitance Ca which correspond to a random charging energy Ec = e^C^^ /2. 

As is well known [1], the QPM (1) is equivalent to the boson Hubbard model (BHM) 
in the limit of large particle numbers per junction. The BHM describes soft-core bosons 
hopping on a lattice [6] and is defined as 

H = iT^riiUijUi - /^E^i - I T^iblbi + H.C.). (2) 

^ ij i ^ {ij) 

Here, b\ (hi) is the creation (annihilation) operator for bosons and rii — blbi is the number 
operator. By writing the field hi in terms of its amplitude and phase and by neglecting the 
deviations of the amplitude from its average, we are lead to the QPM (1). An exact mapping 
between the two models has been derived in Ref. [8]. The hopping term is associated with the 
Josephson tunneling {{n)t — > Ej) whereas C/jj — > 4e^Cy ^ describes the Coulomb interactions 
between bosons. The chemical potential in the BHM plays a role analogous to the one of the 
external charge in the QPM {/j, — > ^i). Thus a QPM with random offset charge corresponds 
to a BHM with random on-site energies. 

The disordered BHM has attracted much attention during the last decade [6,9-14]. In 
the pioneering work of Fisher et al. [6] the phase diagram at T = was studied. Without 
disorder, the BHM exhibits two types of phases: a superfiuid phase and a Mott insulating 
phase, with the latter characterized by integer (or commensurate) boson densities, by the 
existence of a gap for the particle-hole excitations, and by zero compressibility. The phase 
structure is determined by the two competing terms of the Hamiltonian: the charging energy 
leads to a charge localization in the array, while the Josephson energy induces a phase 
coherence giving rise to the overall array superfiuidity. When disorder is present, a third, 
intermediate, phase occurs: the Bose glass (BG). This phase has an infinite superfiuid 
susceptibility, but no gap and finite compressibility. It was early realized that considering 
uniform probability distributions there is no BG in the MF [6]. However, it has been shown 
that introducing a not uniform triangular probability distribution of offset charges the BG 
phase appears in the MF phase diagram [9] . A direct Mott insulator to superfiuid transition 
without an intervening BG phase for weak disorder has been recently investigated in two 
dimensions [10,11]. Very recently, the phase diagram of the two-dimensional disordered 
BHM has been studied at T = by means of Monte Carlo simulations [14] evidencing the 
existence of a Bose glass to superfiuid transition in the strong-disorder regime. 

In the present paper we shall consider uniform disorder distributions: our MF approach 
distinguishes the phases with order parameter (cos ip\) equal to zero (insulating region) or 
different from zero (superconducting region), but it cannot capture the BG. Therefore, the 
full phase diagram is much richer than the MF one. It is obvious that for d = 1 (e.g. JJ 
chains [15]) the MF fails to provide reliable information. However, the MF allows us to take 
a first step towards the understanding of the properties of the disordered JJA and BHM at 
finite temperature. Furthermore, conventional wisdom suggests that for large dimensions 
the MF approach provides the correct phase diagram of the system. 
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The plan of the paper is the following: in Section II, we outline the MF theory for the 
pure quantum JJA and then we compute the free energy averaged over the disorder. In this 
way we are able to get a general formula for the phase boundary line at finite temperature. 
Section III is devoted to the study of the effects of random offset charges with diagonal and 
nearest-neighbor capacitance matrices. In Section IV we consider random self-capacitances. 
To check our results, we study in Appendix A the infinite-range hopping limit and provide 
in Appendix B an alternative and more intuitive MF approach for JJA's in the presence of 
capacitive disorder. Finally, Sec. V is devoted to concluding remarks. 



II. AVERAGE OVER THE DISORDER IN MEAN-FIELD THEORY 



In a functional approach which makes use of the Hubbard-Stratonovich transformation, 
the partition function of quantum JJA's may be written as [5] 



(3) 



where (3 = l/ksT and 7ij = 1 if i, j are nearest neighbors and equals zero otherwise (i.e., the 
hopping term just between nearest neighbors). If 7ij = 1 for all pairs i, j on the lattice, we 
are led to the infinite-range hopping limit which provides a remarkable example of exactly 
solvable MF theory [6] . In the following we shall treat exphcitly also the infinite-range case. 
In Eq. (3), Sf.ff is the effective action for the auxihary Hubbard-Stratonovich field ■^i^ 

SEff M = - In j/ n exp E + ^ E(?i| - V'ie'^' - Vre"^^') 

(4) 

The field V'i may be regarded as the order parameter for the insulator-superconductor 
phase transition because it turns out to be proportional to (e"^'). Since the phase transition 
is second-order [16], close to the onset of superconductivity the order parameter is small. 
One may then expand the effective action up to the second order in ■^i, getting 



(5) 



F[ip] is the Ginzburg-Landau free energy, which - to the second order in ipi - is given by 



ij 







(6) 



where the phase correlator is given by 

1 S^Sf 



/?2 Si:i{T)Sijj{T') 



(7) 



A straightforward, but lengthy, calculation gives [5] 
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Y, , e-Eij2e^/3CiTi(ni+|i)(nj + ^)-Ei4e2CTiK+|i)(r-r') 



(8) 



Here, rii assumes all integer values and J2[ni] ^ over all the configurations. If one 
introduces the Fourier transforms in the imaginary time and in the space 



1 



(9) 

(10) 



(with LOfj, Bose-Matsubara frequencies and k vectors of the reciprocal lattice) , the Ginzburg- 
Landau free energy (6), reads [7,17] 



2 Gk-k'(^A- 
E]^^ '^•^'^ N- 



(11) 



where 



^ij jy X] Tk 



(12) 



7j7^ is the inverse of the Fourier transform of the Josephson coupling strength 7ij; since 
7ij = 1 if i, j are nearest neighbors and zero otherwise, we have = (Sp e"*'''^)"^, where 
p is a vector connecting two nearest-neighbors sites. Expanding in k one gets 7^^ = l/z+- • •, 
where z is the coordination number. Substituting in Eq. (11) and keeping only the lowest- 
order terms in uj^, k and 1/z, the MF Ginzburg-Landau free energy reads 



— E 



Ejz 



(13) 



In Eq. (13), Gq is 



<^o = ^E<^r(a;/^ = 0'^ = ^c). 



(14) 



As evidenced in Ref. [5] the MF theory approximation amounts to neglect all higher order 
terms in Eq. (13). 

For a given realization of the disorder, the Ginzburg-Landau free energy (i.e., the free 
energy near the transition) is given by Eq. (13). The average of the free energy over all 
the possible reahzation of the disorder allows for evaluation of the effect of a random charge 
frustration {q\} or a random diagonal charging energy Ua — 4e^Cjj"^: one has then 



= / d{X}P{{X})F[ 



(15) 



where P{{X}) is a given probability distribution and d{X}P{{X}) = Y[idqiP{qi) when 
one considers random offset charges or d{X}P{{X}) = YlidUnP{Uii) for random charg- 
ing energies. The random variables on different sites are talcen to be independent. The 
phase boundary line between the insulating and the superconducting phase is determined 
by requiring that F — 0, which in turn leads to 

1 = z^Go. (16) 

In Appendix A we will discuss how Eq. (16) is modified in the infinite-range hopping limit. 
A comparison of the MF discussed in this Section with an alternative MF approach in the 
presence of disorder is given in Appendix B. 



III. RANDOM OFFSET CHARGES 

In the following, we shall consider three different random offset charges probability 
distribution with mean q and width a. That is, a Gaussian distribution P{qi) = const 
X • e~*^^'~^^ /^'^ , a uniform distribution P{qi) = const between q — a, and q + a and a sum 
of 5-like distributions P{qi) — ^nPn^iqi — ne), with J2nPn — 1- For a diagonal capacitance 
matrix, Eq. (16) leads to 

- = / dqP{q)g{q,y), (17) 
a J 

with a = zEj/4:Ec, y = kBTc/ Ec, and 

V- -^(»+(//2e)^ 1 

We observe that the lobes of JJA are invariant under q/2e — > q/2e + l and symmetric around 
q/2e — n + 1/2, where n is an integer: 

g{q + 2ne, y) = g{q, y); g{n + 1/2 + q/2e, y) = g{n + 1/2 - q/2e, y). (19) 

If one considers the infinite-range hopping hmit, one still gets Eq. (17) provided that 
a = J/4:E^. 

The results obtained from Eq. (17) with a Gaussian distribution are displayed in Fig. 
1. One observes that when q/2e = 0, increasing a favors the superconducting phase while, 
when q/2e = 1/2, increasing a leads to the increase of the insulating phase. For large a 
(i.e., (7 > e), the phase boundary hue is the same for all the values of q (in Fig. 1 the 
large-(7 behavior is represented by the bold line). This is expected since, when a is large, 
the average free energy F does not depend any longer on q. 

A useful representation of the phase diagram is obtained by plotting, at fixed tem- 
perature, the phase boundary line on the plane q-a. Without disorder one observes the 
well-known lobe structure [1]. In the presence of weak disorder and at T = 0, the lobes 
shrink, evidencing a decrease of the insulating phase: for a Gaussian (or unbounded) dis- 
tribution the insulating phase completely disappears even for an arbitrarily weak disorder 
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[6]. In Fig. 2 we exhibit the phase boundary hne on the plane q — a at finite temperature 
for the Gaussian and uniform distributions. Of course, even at finite temperature, when 
the disorder increases the lobes flatten and the same lobe structure is obtained from both 
distributions. At T = we recover the result of Ref. [6]. This can be easily seen if one 
observes that, at low temperatures, for 1^1 < e, one has [7] 

^(?,y-o) = ^-i^. (20) 

Without disorder (a = 0), Eq. (17) simply gives a = 1 — 4(g/2e)^. With the Gaussian 
distribution, since g has a pole in the half-integer value of the Cooper charge, the integral 
in Eq. (17) diverges, and a — > 0; i.e., the lobes disappear for every value of a. As evidenced 
in Fig. 3, for a uniform distribution, when a > e, then o; — > 0; when o" < e, a — > only for 
e — a<q<e + a in agreement with [6] . 
We now consider the case 

^(?i)=E^^-%i-^e), (21) 

n 

with Y^nPn = 1; corresponding to a random distribution of charges which are integer mul- 
tiples of e. Actually, this is the most realistic situation for a random distribution. Indeed, 
the probability distributions employed before should be viewed as fictitious continuous dis- 
tributions; i.e., the properties of the overall distribution of charges (mean value and width) 
can be well approximated with a continuous distribution P{q). Inserting the probability 
distribution (21) in (18) we have 

dqP{q)giq,y) = / dq'^pj{q - ne)g{q,y) = ^Png{ne,y) + ^Png{ne,y) 

n odd even 

where J2odd (J2even) IS a sum restricted to odd (even) integer. Recalling Eq. (19), we have 
g{2ne,y) = g{0,y) and g{{2n+ l)e,y) = g{e,y), and we find 

-^Pog{o,y)+Peg{e,y) (22) 

a 

where po = J2evenPn {Pe = J2oddPn) IS the probability that the offset charge q is an even 
(odd) integer multiple of e. In Fig. 4 we plot the phase boundary line (22) for po = Pe = V^- 
We are able to show now that applying the MF approximation described in Appendix B 
with the probability distribution (21) it is possible to find exactly Eq. (22). The eigenvalue 
equation for the Hamiltonian (Bl), Hiipn = En'4^ni can be recast in the standard form of the 
Mathieu equation with the transformation ipnifi) — ^~^^'^'Pn{^i)- We find 

__ +^___ eosv.jp. = 



where A. = En/Ec and v = —zEj{cos(p) /2Ec- The periodic boundary condition ^jjni^i) = 
ipni^i + 27r) gives Pni^i) = PniVi + 27r) for qi = 2ne and Pn{^i) = -Pni^i + 2^) for qi = 
2(n + l)e (n integer). For small values of v, the periodic and antiperiodic solutions of the 
Mathieu equation can be calculated analytically, giving exactly Eq. (22). 
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A. NONDIAGONAL CAPACITANCE MATRICES 



With nondiagonal capacitance matrices, the phase diagram without disorder becomes 
richer [1]. For concreteness, we shall consider on-site and a weaker nearest- neighbor (NN) 
interaction; i.e., the inverse capacitance matrix is restricted to diagonal and NN terms. If 

one defines 9 as the ratio between NN and diagonal terms, one should restrict only to 2;^ < 1 
in order to insure the invertibility of the capacitance matrix [18]. Without disorder, at T = 
an insulating lobe around q = e appears: the width of this lobe is z6/{l + z6). Putting 
W ^ 1 + ze, Eq. (18) for \q/2e\ < 1/2W gives g{q,y ^ 0) = 1/[1 - AW^{q/2ef]; for 
1/2W <q/2e<l- {1/2W) it becomes [7] 



9{q,y^ 0) = - 



1 



+ 



{2Wj- - l){2Wj- - 3) {2W{f^ - 1) + l){2W{f^ - 1) + 3) 



(23) 



In presence of disorder, Eq. (17) for a uniform distribution gives a = for {\/2W)—a < q < 
{1/2W) + (7 and 1 - {1/2W) -a<q<l- {1/2W) + a. Thus, the lobe width decreases as 
{z6 — 2aW)/W. One sees that for a = zd/2W the insulating lobe around q — e disappears. 
This phenomenon is evidenced in Fig. 5. 



IV. RANDOM SELF-CAPACITANCES 

In the following we limit our analysis only to JJA's with random self-capacitance Ca 
and uniform charge frustration q. We shall consider, in fact, a random diagonal charging 
energy U\\, which we assume to be independently distributed according to the probability 
distribution P{Uu) oc e^*^^"^^") Z^"" , where the diagonal electrostatic contribution to the 
energy Ua needs to be positive. By averaging the free energy (15), the equation for the 
phase boundary becomes 

1 P(U) 

-= / dU^g{U,y) (24) 
a Jo U 

where now a — zEj/AUq and the function g{U, y) is given by 

The results of Eq. (24) are summarized in Figs. 6 and 7: when a is small, the supercon- 
ducting phase increases in comparison to the nonrandom case: this is due to the factor 1/U 
in Eq. (24), which makes larger the contribution of junctions with charging energies less 
than Uq. The increase of the superconducting phase is thus due to a decrease of the effective 
value of Ec- This behavior occurs until a reaches a critical value (depending on the charge 
frustration and on the temperature), of order Uq: at this value of a the insulating region 
starts to increase. This is due to the asymmetry of the distribution, which has its peak in 
Uq, but only for positive values. This phenomenon is present also if one considers different 
distributions. An interesting observation is that, when q/2e = 1/2 (maximum frustration 
induced by the external offset charges), the randomness does not modify considerably the 
phase diagram. This should be compared with the nonfrustrated case {q/2e — 0), where 
randomness sensibly affects the phase diagram. 
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V. CONCLUSIONS 



We obtained the phase diagram at finite temperature of JJA's with capacitive disorder 
(i.e., random offset charges and/or random charging energies) by using a MF approximation. 
For a random distribution of offset charges with mean q and variance cr, one has that for 
cr ^ e, the phase boundary fine coincides for any value of q and the lobe structures on 
the plane q ^ a disappear (a is the ratio between the Joscphson and charging energies). 
At T = the result of Ref. [6] are retrieved. If one considers in the BHM also a nearest- 
neighbor interaction, the insulating lobe around q — e, which arises in absence of disorder, 
is destroyed even for small values of a. For arrays with random charging energies, when the 
variance of the probability distribution is smaller than a critical value, the superconducting 
phase increases with respect to the situation in which all self-capacitances are equal. To 
check our results, we have considered also the infinite-range hopping limit of the QPM. 

Within the MF approach used here, it is not possible to capture the Bose glass phase 
and therefore the full physical picture can be much richer than the one extracted from 
MF. However, if in low-dimensional systems (e.g., in Josephson junction chains) the MF is 
expected to fail [15], we envisage that our results provide qualitatively correct predictions 
in large dimensional arrays. 
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APPENDIX A: INFINITE-RANGE HOPPING LIMIT 

We consider the Hamiltonian (1) in the limit of infinite-range hopping, i.e., 7ij = 1 if 

i 7^ j and 7ij = if i = j. A meaningful thermodynamic limit (N — > oo) is ensured by the 
scaling of the Josephson term Ej = J/N. This model has been studied in Ref. [6] for the 
disordered BHM at T = 0. Here we study the finite temperature case within the approach 
provided in this paper. From 7ij = Xlk 7ke''''*^'~-'-', it follows that 7k = iV^k,o ~ 1- Therefore 
the free energy (11) is 



1 



2NS- 



kk' 



G 



k-k' 



N 



'0k'(<^/i)- 



(Al) 



Averaging over the disorder, requiring F = 0, for A?" — > oo one has 



(A2) 



A comparison with Eq. (16) makes evident that the results of Sees. HI and IV apply 
also in the infinite-range hopping model, provided that zEj — > J. 
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APPENDIX B: AN ALTERNATIVE MEAN-FIELD APPROACH 



The effect of the quantum fluctuations is generally underestimated in MF theory: we 
expect that a MF theory approach could provide qualitative information on the phase dia- 
gram in three dimensions, but not in lower dimensions. In fact, a similar situation arises in 
quantum spin glasses where MF treatments are able to establish the existence of a spin glass 
transition [19]. Even if powerful and elegant, the mean-field approximation in the functional 
approach should be regarded only as a first step in understanding the role played by the 
disorder in the mean- field average which leads to Eq. (16). 

The aim of this appendix is to make more explicit the relationship between the average 
over the disorder and the quantum statistical average by resorting to a molecular field theory 
approach. In fact, we shall present an alternative and more intuitive mean-field approach, 
which may be easily applied only to JJA's with diagonal capacitance matrices. 

For this purpose, let us consider a JJA's with a fixed charging energy Eq and a random 
distribution of offset charges -P(gi). The charging term in the QPM Hamiltonian (1) is 
then diagonal and the Josephson term couples different sites. Mean-field theory consists in 
replacing the Josephson coupling of the phase on a given island i to its neighbors by an 
average coupling: EjJ2{ii) cos(93i — ipj) = zEj{cosip) X^iCos^^i. In this way the Hamiltonian 
becomes a sum of single site Hamiltonians, H = J2i Hi, where 



The single-site Hamiltonian (Bl) depends on the random charge frustration qi: thus, 
its eigenf unctions and eigenvalues depend on qi. Therefore, one has to impose the self- 
consistency condition with a double average, the quantum one and the average over the 
disorder: 



where the ipn arc the cigcnfunctions of the single-site Hamiltonian Hi . The phase boundary 
line is obtained from Eq. (B2) by requiring (cosy?) to be small and by keeping only terms 
proportional to it (wc recall that the transition is second order). 

We conclude this appendix by stressing that the application of mean- field theory in the 
presence of disorder corresponds to introducing an order parameter, which is averaged also 
over the disorder: the self-consistency condition then gives the correct mean-field phase 
boundary line. In Sec. HI, Eqs. (B2) and (16) are compared for the 5-like probability 
distribution, showing that they give the same results. 



Hi = -4:EcTr^ - SiEc 



2e difi 



zE J {cos (p) cos </?i. 



(Bl) 




(B2) 
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FIGURES 



T 

c 




FIG. 1. Phase diagram for random offset charges with Gaussian distribution (Tc is in units of 
ks/Ec). The bold hne is for a = e and it represents the case of large variance and it is found for 
all values of the mean q. To the left (right), we plot q = e {q = 0); we use a/2e = 0.1 (dotted line) 
and 0.25 (dashed line). The I and S indicate, respectively, insulating and superconducting phase. 



(a) (b) 




^%.0 0.5 1.0 0.0 0.5 

q/2e q/2e 

FIG. 2. Phase diagram with diagonal capacitance and random offset charges with Gaussian (a) 
and uniform (b) distribution. Top (bottom) of the figures: ksT/Ec = 1(0.1). We plot the cases 
a/2e = 0.1 (solid lines), 0.25 (dotted lines), 0.5 (dashed lines). For large a/2e the phase boundary 
line is flat and it is the same for both distributions. 
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FIG. 3. Phase diagram at T = for random offset charges with uniform distribution and 
short-ranged inverse capacitance matrix. a"/2e is respectively 0.1 (sohd hne), 0.25 (dotted hne) 
and 0.40 (dashed Hne). 
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FIG. 4. Phase diagram for random offset charges with probabihty distribution given by (21) 
and short-ranged inverse capacitance matrix. In the plot Pq = Pe = 1/2. 
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FIG. 5. Phase diagram at T = for random offset charges with uniform distribution and 
short-ranged inverse capacitance matrix. In the plot z9 is equal to 0.1 while a/2e is respectively 
(solid line), 0.015 (dotted line) and 0.03 (dashed line). For this value of z9, the lobe disappears at 
a = 0.045. 




FIG. 6. Phase diagram in the Tc — a plane for random diagonal capacitance with Gaussian 
distribution and uniform offset charge q/2e = (a), 0.5 (b) (Tc is in units of kB/Uo) while cr/2e is 
0.1 (solid line), 1 (dotted line) and 5 (dashed line). 
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FIG. 7. Phase diagram in the q — a. plane for random diagonal capacitance with Gaussian 
distribution and uniform offset charge at ksTj^Q = 0.1. (T/2e = 0.1 (solid line), 1 (dotted line) 
and 5 (dashed line). 
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